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Abstract 

The classical van Royen-Weisskopf formula for the decay width of a meson 
into a lepton-antilepton pair is modified in order to include non-zero quark 
momentum contributions within the meson as well as relativistic effects. Be- 
sides, a phenomenological electromagnetic density for quarks is introduced. 
The meson wave functions are obtained from two different models: a chi- 
ral constituent quark model and a quark potential model including instanton 
effects. The modified van Royen-Weisskopf formula is found to improve sys- 
tematically the results for the widths, giving an overall good description of 
all known decays. 
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I. INTRODUCTION 

In a previous paper |l[ we studied meson strong decays in the framework of non- 
relativistic quark models. In that work, a comparison was made between a potential in- 
cluding only gluon exchanges and a potential based on a chiral constituent quark model 
(xQM). Although originally designed for the nucleon-nucleon (NN) interaction, the xQM 
has been successfully applied to the baryon and meson spectra, besides the meson strong 
decays. This is a stringent test for the model, as there are only a few free parameters, which 
are taken to be the same for the NN interaction and the hadronic properties. Such a unified 
description reveals also the capability of a non-relativistic approach to take care of a wide 
set of phenomena, even in situations where it seems not justified to apply it. 
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The calculation of the spectrum is just an average of the meson wave function; therefore, 
to discriminate between different models it is necessary to rely on more sensitive observables. 
Meson strong decays were already calculated and, as already said, with rather good results. 
But electromagnetic transitions are even more interesting, since the transition operator is 
perfectly known and therefore they constitute a very drastic test of the wave functions. The 
present work is a first step to the study of electromagnetic properties, it deals with the decay 
of a meson into a lepton-antilepton pair. 

Since many years the study of the decay width of a meson into a lepton-antilepton pair 
has been done in the frame of the van Royen-Weisskopf formula |J. This work is based 
on a second order Feynman graph, where the quark- ant iquark pair which forms the meson 
annihilates into a virtual photon which later on produces the lepton-antilepton pair. This 
implies that only neutral mesons with J p = are allowed to decay within this scheme. In 
practice, this concerns 3 Si and 3 Di resonances (strictly speaking a true state can be a mixing 
of both, due to the tensor force, but in this study we discard it so that the corresponding 
states are decoupled). The calculation is realized through the familiar Dirac trace method 
in order to evaluate the squared modulus of the transition amplitude. In the original work 
0, besides a missing factor 3 due to the ignorance of colour degrees of freedom at that 
time, there are mainly two approximations: first of all to consider point-like (bare) charged 
quarks, and secondly to consider the quarks having zero momentum inside the mesons. 

The approximation of point-like charged quarks is not consistent with the assumptions 
of the constituent quark model. Constituent quark masses arise when an initially massless 
quark moves through the QCD vacuum. Due to the strong qq interaction, the quarks get 
dressed with a polarization cloud of qq pairs. Therefore, a finite electromagnetic size appears 
due to the quark-antiquark dressing. 

The second approximation is the responsible for the leptonic decay widths being deter- 
mined by the value of the meson wave function at the origin. This quantity is very sensitive 
to the underlying potential, so that meson leptonic decay widths are really a very impor- 
tant observable to discriminate between different wave functions. Within the original van 
Royen-Weisskopf formalism, 3 Di states, having a null wave function at the origin, are not al- 
lowed to decay. Nevertheless, some decays of such states have been observed experimentally, 
indicating limitations to the applicability of the van Royen-Weisskopf formula. 

There have been other studies of the leptonic decay widths of mesons. In Ref. the 
leptonic widths of light vector mesons were calculated in a modified MIT bag model under the 
assumption that the total quark-antiquark momentum is zero in the center-of-mass system 
of the meson. Ref. |4| used a similar formalism but based on the dominance of a scalar- 
vector harmonic confining interaction for the mesons. Both calculations neglected residual 
interactions like one-gluon exchange. Extensions of the van Royen-Weisskopf formalism can 
also be found in the literature. The Bonn group || obtained a quasi-relativistic formula 
by considering p/m terms in the decays, i.e., without neglecting the small components of 
the Dirac spinors of quarks. The modification to the van Royen-Weisskopf results for S- 
waves ranges between 5 and 15 % depending on the state. A non-zero, although very small, 
probability for D-w&ve meson decays was obtained. Ref. |J also considered relativistic 
corrections to the van Royen-Weisskopf formula in a covariant formulation of light-front 
dynamics. The difference with respect to non-relativistic calculations can be as large as 60 
% in the charm sector, being smaller (~ 30 %) for bottomonium. 
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Our purpose in this work is to extend the van Royen-Weisskopf formula in order to 
take care of the a priori important aspects mentioned above. First of all, to consider 
the momentum distribution of quarks inside a meson, as well as kinematical relativistic 
effects in the non-relativistic wave function. Secondly, to study the role played by the 
constituent quark electromagnetic size. We will work in the framework of the constituent 
quark model of Ref. J7| assuming that the meson decay proceeds via the qq annihilation 
as in Ref. ||. Therefore, our formalism will be on the same spirit as that of Ref. [[J. In 
order to calculate the corresponding Feynman amplitude we modify our non-relativistic wave 
function including a relativistic spinor. Finally, we want to test the sensitivity of the results 
to the corresponding non-relativistic wave functions. For this purpose we will perform an 
exhaustive study of all seen transitions using the two different quark-antiquark potentials 
previously mentioned. 

The paper is organized as follows. In the next section we calculate the decay width of 
the transition taking into account all the improvements presented above. The main features 
of the quark-antiquark potentials are refreshed in section |T|. The results are presented and 
analyzed in section [IV]. In the last section we summarize the conclusions of our work. Some 
technical aspects as well as the specific form of the electromagnetic size for the quarks are 
detailed in the appendices. 

II. CALCULATION OF THE DECAY WIDTH 
A. Notation 

As the formalism is very involved let us first of all make explicit the notation we will use. 

For the meson: mass M, momentum in the center-of-mass reference frame P = 0, 
total angular momentum J and projection on a fixed z-axis fij, orbital angular momen- 
tum L and projection fii, spin S and projection /is, wave function in coordinate space 
^(r) = R(r)YLfj, L (9,ip), and wave function in momentum space 4>{p) = R p (p)YLfj, L (9,ip) (a 
more correct expression for the arguments of the last spherical harmonics would be (8 P , ip p ), 
however for typographical facility we still use (6,ip) when no confusion arises). 

For the quark: mass m q , charge e q (in units of the elementary positive charge e), mo- 
mentum p, kinetic energy E p = ^Jrn 2 + p 2 , spin wave function \&m) with projection M, and 
helicity wave function \x\) with projection A. 

For the lepton: mass m;, charge — e , momentum q, kinetic energy E q = \Jm 2 + q 2 (do 
not confuse this q index with the index relative to the quark; in principle the context is not 
ambiguous) and spin projection £. 

For the antiparticle (antiquark and antilepton): the mass is the same as the one of the 
particle, the charge is opposite and, since we are in the rest frame, the momentum is opposite 
and the kinetic energy identical. 

In addition, the virtual photon four-momentum p 7 verifies p 2 = M 2 by energy conserva- 
tion, and we express the Clebsch-Gordan coefficients as {iimij2m2\j^m^) . 
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B. Leptonic decay width 



We calculate the decay width in the S matrix formalism. One should remember that 
the Feynman diagram refers to free particles with given momenta; therefore, in a meson, we 
must take care of the fact that the fundamental amplitude for the quarks must be integrated 
with the corresponding amplitude probability <fi(p) to find a quark with momentum p inside 
a meson. In the center-of-mass reference frame, the S matrix at second order takes the form: 



' (2?r) 7 



,/V ~"' ' {2nfE p {2nfE q 
11 

2 2 1 ^^^ J M 2 



£ (L^SnsU^^-M^-M^S^s) d 3 p<f>(p) ^A?A» q , ( I ) 

M 1 M 2 (j. l ij,s 

where Af = [ui(q, ti)l^ v i{~ q, £ 2 )] denotes the leptonic part of the S matrix and A v q = 
[■y g (— p, M 2 )Y u q(p, Mi)] the quark part. In this expression we used the fact that the relative 
momentum p appearing in the meson wave function is chosen to be the same as the quark 
momentum and opposite to the antiquark momentum. 

Let us first compute the quark term. Taking for the Dirac spinors the expressions 



t _ E p + m q t ff-p t 
V °-' j 2m q { - M >E p + m q >- M > 



we obtain 



Ep + mq ( $Af, \ 



In the case v — the result is identically zero. For technical simplicity, we switch from 
cartesian to spherical coordinates, so from now on the latin indices will run to — 1, 0, 1. 
Using the relation (a ■ p)<7i(a • p) = 2pi(a ■ p) — p 2 a,i the quark part of the S matrix can be 
written as, 



Ai ^ E p + m q 



2m q 



(E p + m q ) 2 J —^ — ^ ~^ {Ep + mq) 2 



V \ P P 



(4) 



As seen in this equation, to evaluate the S matrix the expectation value of the Pauli 
matrices between spin eigenstates ($m) is needed. These functions are not eigenstates of 
the Dirac free hamiltonian. Moreover, they are eigenfunctions of o"o only for small values of 
the momentum p (but the contributions of large momenta are suppressed by the exponential 
tail of the wave function). Nevertheless, the calculation of A q can still be done with the $m 
spinors, but being very careful with phases. Alternatively, one can switch to the helicity 
scheme, for which the spinors are eigenstates of the Dirac free hamiltonian. We will work in 
the last scheme, but we checked that both methods give exactly the same results. 
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Given a spin one-half particle with momentum p = (p,9,ip), the relation between spin 
and helicity eigenstates is M: 



\®m) = Y, D -1-m&9,<p)\Xx(0,<p)), 



^M\ = T,(-^ x ~ M Dl%(o,e, v )(xx(e, v )\, 



(5) 



where (0,9, if) are the Wigner matrix elements. For a two-particle state with mo- 

menta p and —p, in order to keep the system invariant under Lorentz transformations, it is 
necessary to introduce an additional phase in the wave function of the second particle ||. 
This phase is (— 1) J_A , being j the particle spin and A its helicity, and therefore the spinors 
take the form 



Ai 

(S-Afcl = E(- 1 ) 1/2+M2 < 2 -M 2 (0^^)(XA 2 (^^)|. 
A 2 



(6) 



Using the Wigner matrix elements and the mean value of the Pauli matrices between 
helicity states given in Table §, and summing over all helicity states, we get the following 
results for $-M a a i$Mi as a function of the spin indices Mi and M 2 , 



(M 1 ,M 2 ) 


(1/2,1/2) (1/2,-1/2) 


-1/2,1/2) (-1/2,-1/2) 




— V2<5j_i 5 i>0 





One can then perform the summation over Mi and M 2 in Eq. ([!]) arriving to 



J2 (lM^M 2 \Sfis)^. M2 a^ Ml = V2(-iy +1 5s,i5 ir 

MiM 2 



MS' 



from where a first selection rule for the leptonic decay: S = 1, is obtained. Eq. 
finally written as 



At = 5 S iV2 



^(-1) J+1 5. 
m„ 



1 



where it has been used that 



1 + 



m q (E p + m q ) 



2 En 



j2PiP*(-iy +1 s 



J, -Ms 



(E p + m q ) 2 E p + m q 



(7) 
can be 

(8) 
(9) 



Let us now concentrate on the part of the S matrix which depends on the quark and the 
antiquark, 



B l = J d 3 P 0(p) £ (Ln L Sn a \Jnj)5 8 ,i\/2 



M-LMS 



(10) 
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The integration of the first term can be done through the relation / d 3 p<p(p) = 
(27r) 3/2 vI/(0), arriving to 

(2n) 3 ^(0)5 L , 5 s ,i5j, 1 5 l ^ J V2(-l) 1 ^^ (11) 

The integration of the second term is a little bit more involved, and after some algebra 
one arrives to 

v / 2v / 2v / 2TTl(27r) 3 / 2 



6tt 

where we have defined 



Taking the two parts together, Eq. fllOf) is finally written as 



■(LOlOllO)/^,!^,!^^-!) 1 -^, (12) 



, p 4 d P R p (p) 

h ~ 1 E p (E p + m q y {U) 



V2 



^(0)h,o ~ j-V2L + l(L010\10)h 

D7T 



(14) 



The 5 function and the Clebsch-Gordan coefficient provide a new selection rule: L = 0, 2. 
The first term between brackets gives rise to the classical van Royen-Weisskopf formula, it 
only concerns to 3 Si resonances. The second term is the modification due to non-zero quark 
momentum distribution inside the meson (through the integration over p in I4) and it also 
contains kinematical relativistic effects (through the E p denominator in I4). It contributes 
to L = ( 3 Si) resonances and L = 2 ( 3 Di) resonances, whose decays are strictly forbidden 
within the original van Royen-Weisskopf formalism. The recipe for taking into account the 
previous improvements is quite simple; one makes the substitution 

V2 



tf(0) -> ^(Q)5 Lfi - j^V2L + T(L010|10)/ 4 = *(0) (15) 

in the traditional van Royen-Weisskopf formula. 
Finally, we have for the S matrix, 

S = ^ 2e ^^^v / 2(2 7 r) 3 / 2 ^(0)5( 4 )(g 1 + q 2 - P) £ ( 16 ) 

The rest of the calculation to obtain the expression of the width is detailed in appendix 
[A[ Up to now, the S matrix has been calculated for a quark- ant iquark system with given 
colour and flavour. Considering that the meson is a colour singlet with a contribution of each 
colour with the same probability, the net colour factor is a multiplication by 3 as compared 
to the usual van Royen-Weisskopf expression. For neutral ordinary mesons composed of two 
flavours / (u and d) with quark charge e q (f), the flavour wave function is written generally 
as Y^fOjfff (the amplitude a/ on the // channel is simply a Clebsch-Gordan coefficient). 
Provided an SU(2) invariance of the potential, the meson wave function is the same for 
the two flavour channels, so that flavour degree of freedom implies only the replacement 
e q — > J2f a f e q(f)- The final result for the width is: 

r = i ^( 1 -4)' /, ( 1 + 4)C«^)'^, (17) 
being a the fine structure constant. 
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C. Electromagnetic size for quarks 



In the previous section we have treated quarks as point-like charge distributions. How- 
ever, the constituent quark is a complicated object surrounded by a quark- ant iquark pair 
cloud. Vogl et al. [[Uj have explicitly shown in a generalized Nambu-Jona-Lasinio model that 
constituent quarks are dressed by a mesonic (qq) polarization cloud. For a given external 
field the qq modes, which carry the quantum numbers of the probe, will determine the size 
of the constituent quark. Hence, the electromagnetic properties of constituent quarks in the 
low-energy regime are largely determined by the screening effects induced by qq polariza- 
tion cloud with p meson quantum numbers, and the old vector-dominance principle emerges 
naturally. 

Based on these results several authors have studied hadronic electromagnetic properties 
assuming that the electromagnetic size of the constituent quark can be described by a 
monopole form factor 

F(q 2 ) = * 2 - (is) 

1 + q A m l 



Magnetic form factors [[LI], magnetic moments [12|], electromagnetic properties of the 



A (1232) [|H|, pion form factors Jl4| and radiative up form factors [15], among others, have 



been calculated in different models. The results show the relevance of the quark electromag- 
netic size. 

The quark electromagnetic size is then related to the corresponding vector-meson prop- 
agator and it has a dependence on the inverse of the squared vector-meson mass. In this 
way, the quark electromagnetic size is flavour- dependent, through the flavour dependence 
of the vector-meson mass. Instead of using the preceding expresion, we simulate the quark 
electromagnetic size by a density p q (r). These quarks will be called from now on dressed 
quarks. Let us denote by a^(f) the creation operator of a bare quark at position r; in the 
same way A^(r) stands for the creation operator of a dressed quark at position r. It seems 
natural to define: 

A\r) = f d 3 r'p q (f-f')a\f'). (19) 



If the density function p q (r) = 5^ 3 '(r) , we have A\r) = a)(r) and we recover bare 
quarks. Strictly speaking, the above equation is incorrect and misleading, in particular 
commutation relations among A'(f) operators are no longer satisfied. Nevertheless, we still 
apply this equality keeping in mind that it makes sense only when the p q (r — f') distribution 
is a sharply peaked function. This is the price to pay to get a phenomenological very simple 
modification to a finite size for the quarks. 

Let us calculate now the Fourier transform of the densities and of the creation operators. 
Defining 



re 



ip.r 



fit (p) = jA^ J d 3 ra\r)e^ 

P~ q (p) = j^f 2 J d 3 rp q (r)e^, (20) 
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from the properties of the Fourier transformations it is immediate to see that: 

A\p) = (2nf 2 p q (p)~a\p). (21) 

We define (forgetting the obvious complications due to other degrees of freedom and 
the corresponding couplings) the creation operators of a meson as a bare or dressed quark- 
antiquark pair creation, respectively: 



S = J d 3 p<p(p) d\p)b\-p), 

D ] = J d 3 p<p(p)A\p)B\-p), (22) 

where W and are the corresponding creation operators for bare and dressed antiquarks, 
respectively. 

Operating, we get 



where we have defined 



= {271 f J d 3 p<p(p) p q {p)p q {-p)~a\0{-p) 

= (2ix) 3 j d 3 pmP~Ma j (p)b j (-P), (23) 

Pqqip) = Pq(p)Pq(-p). (24) 



Note that, since the antiquark has the same mass as the quark and the densities depend 
only on p 2 , it is natural to take p q (p) = Pq{—p), so that: 

P~M=p q (p)\ (25) 

which implies 

Pqqip) = P~qq(-P)- ( 26 ) 

When operating on the vacuum state |0), gives the spinors u q and v q of Eq. @. 
In the previous section, we arrive to Eq. fllPP, in which there are two terms. The only 
dependence of the first term in p is in the wave function; with dressed quarks we would have 

?3, 



(27r) 3 y d'pmp^M = w y d^mpU-p) 

= (2n) 3 J d 3 p<P(p) j^j- 2 J d*re-^p qq -(r) 

= (27r) 3 Jd 3 r^(r)p qq (f), (27) 

while the second term would simply include the quark density. 

Therefore, when introducing dressed quarks, it is necessary to substitute: 

*(0) -> (2tt) 3 / 2 J d 3 rp qq (r)^(f), (28) 
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III. MESON WAVE FUNCTIONS 



Once the mechanism to calculate the decay widths has been described, the only remain- 
ing point is the meson wave function. The sensitivity of the results on the meson wave 
function will be analyzed by comparing two different potential models which seem to us 
very representative of the high degree of sophistication reached in meson spectroscopy. The 
first one is the chiral constituent quark model, xQM. The central feature of the model is the 
concept of constituent quark mass, which appears as a consequence of spontaneous chiral 
symmetry breaking. The potential, derived and discussed elsewhere |],|7[], includes a con- 
fining interaction, a perturbative one-gluon exchange and Goldstone-boson exchanges. It 
can be easily extended to the study of mesons in which at least one of the components is 
a c or b quark (or antiquark). In this case, chiral symmetry is explicitly broken and the 
potential only includes confinement and one-gluon exchange. The parameters used are the 
same as in Ref. |TJ], with the addition of the charm and bottom masses: 1845 and 5250 
MeV, respectively. This potential model provides reasonable results for the baryon-baryon 
interaction ]7), baryon spectra meson spectra ]17| and meson strong decays As a 
consequence, the calculation of meson leptonic decays constitutes a severe test of the model 
as it supposes a parameter-free calculation. 

The second potential, called DNR, has been developed around gluon exchanges JTB) and 
relies very closely to the fundamental QCD. The short-range part sticks to the experimental 
a s (q 2 ) coupling constant while the long-range confinement is essentially linear with a string 
tension roughly equal to the experimental one. Moreover, it includes instanton induced 
effects, which play a crucial role in explaining the isoscalar spectra, and hyperfine terms as 
a sum of two gaussian functions. The quarks have also been dressed with a gluonic cloud 
(which has nothing to do with the electromagnetic cloud considered above), so that the bare 
potential is convoluted with a gluonic density to give the final potential. The free parameters 
have been fitted in all mesonic sectors and the corresponding spectra are quite good. 

Both models rely on a non-relativistic expression of the kinetic energy operator and 
they need to solve the Schrodinger equation (for example using a Numerov algorithm). To 
compare their respective merits we have plotted in Fig. [TJ the spectra for the resonances 
appearing in this study. The quality is more or less the same, some resonances being better 
reproduced with xQM and others with DNR. In particular, except for the T, which raises a 
problem to xQM, both potentials are able to reproduce the heavy sector (charmonium and 
bottomonium) rather nicely, even for very excited levels (up to 6 3 Si in bottomonium). In 
the light sector (ordinary and strange) the situation is a little less favorable; in general xQM 
reproduces better first radial excited states, while DNR looks better for orbital excitations. 

Although both potentials give spectra of roughly comparable quality, they can provide 
rather different wave functions. As an illustration of this fact, we draw in Fig. || the wave 
function R(r) of the p resonance. xQM wave function shows a hole near the origin (due to 
a stronger hyperfine term) which is absent in DNR wave function. As a consequence there 
is a factor « 3 difference for the value of the wave function at the origin (and thus a factor 
around 10 for the width). We thus see the fundamental importance to rely on observables 
other than energies to test the quality of a model. 
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IV. RESULTS AND DISCUSSION 



The previous formalism is applied to all experimentally known decays of a meson into a 
lepton-antilepton pair. Some other transitions have been seen experimentally but without 
precise data; we present also our results in this case as a matter of prediction. 

In order to isolate the dynamical effects from the kinematical ones, we employ the ex- 
perimental meson masses for the phase space factor. Let us emphasize that the calculation 
of the meson leptonic decay widths in the bare quark case is parameter-free, all the model 
parameters being fixed to other observables (meson spectra for DNR, and NN interaction 
and meson spectra for %QM). The results dealing with the dressed quarks need to determine 
the size a of the electromagnetic distribution. As detailed in appendix ||, we use a gaussian 
charge density and the flavour dependence of the size is parametrized by 

a f = a e m f /mo , (30) 

where m/ is the constituent quark mass of the corresponding flavour /, and the parameters 
o"o and mo are fitted on the known transitions. The obtained values are (<7o, mo) = (0.6 
GeV, 1.25 GeV) and (0.7 GeV, 2.11 GeV) for DNR and xQM, respectively. It is worth to 
notice that the electromagnetic sizes coming from this fit are in good agreement with the 
prediction of vector-meson dominance model. 

The results for leptonic decay widths are presented in Table [H[ We compare the value 
obtained for the xQM an d DNR meson wave functions in three different cases. The first two 
columns report the results derived from the original van Royen-Weisskopf formula. Columns 
three and four show the results when the modified van Royen-Weisskopf formula of Eq. ([15]) 
is used. Finally, columns five and six contain the results when, on top of the modified van 
Royen-Weisskopf formula, the dressing for the quarks is used. 

Let us comment first the values due to a pure van Royen-Weisskopf formula ( VR) . In this 
case, we discard the decay of Z D\ resonances that is strictly forbidden. The DNR potential 
systematically gives a too large wave function at the origin; as a consequence the theoretical 
width can exceed the experimental one by a factor 3, but in the bottomonium sector the 
values are much closer to the experimental data. Moreover, the hierarchy in the relative 
importance is respected, indicating VR contains already a part of truth. More puzzling is 
the case of xQM potential. An important difference can be established between the decay 
results for light mesons (p, u and 0) which are smaller than the experimental data and heavy 
mesons (ip and T) which are in general larger than the experimental ones. This is obviously 
due to the hole in the wave function near the origin (see Fig. 0). Although all of them 
have 5=1, and therefore the spin-spin contact term included in the one-gluon exchange is 
repulsive, for light mesons this term is strong, as it depends on the inverse product of quark 
masses, and makes the wave function to have its maximum at r w 0.5 fm. As a consequence, 
the van Royen-Weisskopf value for these decay widths is smaller than experiment. On the 
other hand, for heavy mesons this repulsion is much less important, and the wave function 
takes its maximum at (nearly) the origin. The van Royen-Weisskopf value for decays of 
these mesons is higher than experimental ones, similarly to the results obtained previously 
in the literature. For the heavy sector xQM gives poorer results than DNR. 

For heavy quarkonia, in both models, our extension of the formalism to include non- 
zero momentum for the quarks and kinematical relativistic effects gives always a negative 
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contribution. The difference with non-relativistic results is between 10 and 30 %, bigger 
than in Ref. || but much smaller than in Ref. ||. The main effect is thus to decrease the 
theoretical value, making it closer to experiment. The same conclusion applies for light 
mesons with DNR wave functions. The case of light mesons with xQM is more interesting. 
The modification is positive, so the transition is enhanced making it closer to experiment 
again. Thus, whatever the importance of the wave function at the origin is, the modification 
presented above always improves the theoretical results. This indicates again that there is 
certainly a part of truth in this explanation. In addition, we get non-zero widths for D-w&ve 
mesons, like ^(3770) and ^(4160), although the calculated widths are much smaller than 
the experimental ones, as also observed in Ref. 0. 

The most significative improvement is produced when the electromagnetic size of quarks 
is considered. As the term included in the van Royen-Weisskopf formula is the most relevant 
to obtain the width, the quark finite size replaces the value of the wave function at the origin 
(which is strongly dependent on the structure of the potential at very small distances) by 
an average on the whole wave function, taking into account the behaviour of the system 
everywhere. The I4 integral is also sensitively changed. For DNR the main effect is to 
still reduce the theoretical values and make them very close to the experimental data. The 
results in this case can be considered as very good, owing to the fact that we have only two 
free parameters at our disposal. For xQM the behaviour depends again on which sector we 
consider. For the heavy mesons, the dressing of the quarks reduces the theoretical value 
and makes it closer to experiment. For light mesons the width is importantly increased, 
being also close to the experimental value. Thus, in any case, the description is improved. 
Although the results of xQM are slightly worse than those of DNR, the conclusion is that 
the dressing is very important and it always improves the situation. 

The dressing does not affect very much the decay of D-w&ve resonances, whose calculated 
widths remain still too low. We think that this point must deserve further considerations 
like the inclusion of tensor forces which are outside the scope of this work. 

Concerning the branching ratios between the different lepton channels (e + e~, p + p~ and 
t + t~) our formalism always gives bigger widths for lighter leptons, while the experimental 
results are in almost all cases the opposite. For the moment, we have no explanation 
concerning this point. 

One remark should be done about the ratio between p and u widths. If both mesons 
were degenerated and had the same radial wave function, this ratio would take the value 
of 9 due to flavour content. Since p is experimentally lighter than u, there exists already 
a small deviation due to phase space factor. One would expect also that its wave function 
would be slightly higher at the origin and the ratio would be closer to the experimental 
value, which is approximately 11. In DNR the wave functions are absolutely identical, so 
that this dynamical effect does not exist. With xQM the situation is different because the 
potential is isospin dependent; nevertheless, we find m p > m w and therefore v I / (0)p < ^(0)^; 
as a result, the calculated value for the ratio is even smaller than 9. 

With respect to the transitions for which precise data are not known, both models give 
results which are relatively close to each other, and which, in any case, respect the hierarchy 
of importance. We can thus have some confidence on them. 

Finally, we would like to show how the quark electromagnetic size we have used works 
in other electromagnetic observables. In Fig. we present the pion charge form factor 
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calculated with the xQM wave function used in this work. The dotted line represents the 
calculation done with bare quarks whereas the dashed line corresponds to the dressed quarks 
calculation. As can be seen, the inclusion of the quark electromagnetic size improves the 
agreement with the experimental data. Although it is not our aim to make any definitive 
statement about this problem we can introduce the effect of boosting on the form factor 
(which has been calculated using rest frame wave functions) using the prescription of Stanley 
and Robson JTP| to take into account the kinematical meson recoil. The result of this boosting 
is shown by the solid line, obtaining a reasonable description of the data. It is worth to 
notice that the calculation has been done in the one-body approximation framework for the 
electromagnetic current, without taking into account the two-body currents necessary for 
fulfilling gauge invariance. However, it is clear that a reasonable agreement with the data 
can be achieved by assuming an effective one body current between constituent quarks with 
non-trivial electromagnetic structure. This conclusion agrees with the calculation of Ref. 
n|, altough in a different framework. 



V. SUMMARY 

In this paper we investigate whether the van Royen-Weisskopf formula for meson lep- 
tonic decays can be improved. We first include the possibility of non-zero momentum for 
the quark- ant iquark pair within the meson as well as kinematical relativistic effects. We also 
introduce phenomenologically an electromagnetic density for dressing the quarks, that we 
choose of gaussian form. We make our study on all experimentally known transitions and use 
wave functions coming from two different non-relativistic quark model potentials describing 
correctly the meson spectra. The first one, xQM, is based on gluon and Goldstone-boson 
exchanges, while DNR contains gluon exchange as well as instanton induced effects. Al- 
though the wave functions coming from both models can appreciably differ, we showed that 
the proposed modifications always go in the right direction. 

The inclusion of non-zero momentum for quarks inside the meson gives an effect of order 
10%-30%, while the inclusion of a quark electromagnetic size improves a lot the situation. 
With two adjustable parameters, our theory is able to reproduce nicely the experimental 
data for transitions coming from 5-wave resonances. Our proposed modifications give non- 
zero width for .D-wave mesons, but seem to be not enough to obtain their correct order 
of magnitude. The fact that they systematically improve the description shows that we 
certainly introduce basically good ingredients. 
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APPENDIX A: LEPTONIC PART AND FINAL CALCULATION OF THE 

DECAY WIDTH 



We start from Eq. (0), making explicit the leptonic part 
v / 2(2vr) 3 / 2 ^(0)5( 4 )(g 1 + q 2 - F) 



. 2 l m l 

%t e<? (27r) 2 E g M 2 



-Eg + 

2m; 



1 + 



^ 1 ^-6-2E(- 1 )H^ $ - 



(E q + mi)\ 



Q-iQ-k 



£ 2 



(£, + mi) 



(Al) 



For the leptonic part, as we are performing a sum over all the final states, the same 
result is obtained when using spin or helicity eigenstates. We take spin states for simplicity. 
Using the Wigner-Eckart theorem, 



the S 1 matrix can be written as 



1+ j 2 E g + m l 

2mi 



(A2) 



(A3) 



(E q + mt) 2 I N 2 S1 2 



2 EC- 1 )^^! 1 *)*-*^ 



Once we have the final expression for S, we have to write it in the form |20 

S = -2tt^ 4 >(£p/-E^)M 



(A4) 



in order to calculate the differential width of a process with one initial particle and two final 
particles as 



dT = 2ir\M\ 2 6^(Y,Pf ~ EpO Ai^ 



92, 



(A5) 



performing an average over initial states and a sum over final states. 
In our case, 



dT = (2n)- 2 5^( qi + q 2 - P)^-\^^ ^ ^ ; 



M 4 



d i q 1 d ?J q 2 



E q + mi 
mi 



(A6) 



3,E 



1 + 



(E q + mi) 2 J N 2 Si 2 



(£"5 + mi 



2 EC- 1 )*^^! 1 *)?-*^ 



where a factor 3 is inferred from the colour part and a factor 1/3 from the 2 J + 1 initial 
states average, being J = 1. 

The next step is to compute the square of the modulus of the leptonic part: 
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E 

4 E 



2 E(- 1 ) fe (^4e 2 |lA:}?-^ 



(E q + mi y'2^2 



(kklW 2 + 



p^E(-i) fc <^ 211 ^ 

1 E(-) fc+i (^46|iA:)(^46|i/)g- fc g; J ^^ J 



(E q + ///,)•' 

4E [ ^ 



A'./ 



(£7, + m,) 2 (E q + mi) k 



(E q + mi)' 



ER fc+ W fc <4<?-^ 



3^ 2 



2E, 



" q 2 + 



(Eg + ITU) 2 (Eg + m Z ) 3 (Eg + 771,) 



8£ 2 



+ 



4m 2 



(£ g + m,) 2 ( J E (? + m / ) 2 ' 
Substituting in Eq. 



rfr = (27r)- 2 5( 4 )( gi + g 2 -P) 



9 4 2 

M 4 



4 1 + 



2EI 

mf 



\y(0)\ 2 d i q 1 d i q 2 . 



(A7) 



One of the two integrals can be done using three of the four Dirac deltas, from the other 
one, we obtain a factor 4ir from the angular part, 



167rm 2 e 4 e 2 r q 2 dq 



An 2 M 4 J q 2 + m 

4 mfe 4 e 2 1 
7T M 4 



<5 M — 2A/g 2 + m 2 1 + 2 



t q 2 + m 2 ^ 
m 2 



Mo)\' 



2\ 1-4 



q 2 dq 
q 2 + 



^ \ Q ~ \ ~~T~ ~~ 777,2 )(■'■" : - 



g 2 + m 2 



M 2 



2 mfe 4 e 2 



1 ^-m 2 



7T M 4 



M 2 



1 + 



2\ 1-4 



M 2 
2mf 



M 2 



mt 



M|*(o)p 



(A8) 



and finally, 



r 



16ira 2 e 2 
M 2 



1-4 



n? 



2\ 1/2 
2 \ 



M 2 



'l + 2^)l*Wl S 



(A9) 



which is Eq. flI7|) without the flavour modification explained in the text. This procedure is 
of course more involved and much less elegant than the trace method, but this is the price to 
pay for using a non-relativistic wave function coupled to a good angular momentum inside 
a formalism that is covariant at the beginning. 
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APPENDIX B: GAUSSIAN ELECTROMAGNETIC DENSITY FOR QUARKS 



The density function for quarks must be, in limit of zero size, a Dirac delta. We can 
choose a gaussian function: 



Pq (rl = Jj- 2 e-^. (Bl) 



7T 



In the limit a q — > oo the size of the quark tends to zero and p q (r) — > 5^(r). The Fourier 
transform of the density is 

Therefore, the quark-antiquark pair density is given by 

P~M = 7^i e %? ' ( B3 ) 



1 <tZ, 



where 



= , (B4) 



^ = jMi . (B5) 

V^ 2 + °\ 
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FIG. 1. Comparison of the DNR and xQM spectra for J p = 1~ resonances discussed in the 
text. The horizontal lines indicate the experimental values of the masses, being rectangles in the 
case of an experimental uncertainty bigger than 20 MeV, the open circles show the values obtained 
with DNR potential and the dark circles those from xQM potential. The p(nn system, 1=1), 
u>(nn system, 1=0), (f)(ss system) families have a common scale whereas ip(cc system) and T(bb 
system) families have their own one. The spectroscopic quantum numbers are indicated close to 
the experimental mass values. 
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FIG. 2. Radial wave functions R(r) of the p meson obtained with DNR (solid line) and xQM 
(dashed line) potentials. 
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FIG. 3. Pion charge form factor calculated with the xQM wave function. The dotted line cor- 
responds to the bare quark calculation. Dashed line corresponds to the dressed quark calculation. 
Solid line includes the effect of boosting. 
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TABLES 



TABLE I. Mean value (xAi l c7 i|XA 2 ) °f the Pauli matrices between helicity states. 
Ai A2 o\ <to a -1 

1/2 -1/2 -kte°*± e iv - S [ n - k= ^ 1 

-1/2 1/2 1=^6 e i v _ sin6 1± ^ e 

-1/2 -1/2 s^e*v -cosfl -^e 
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TABLE II. Results for the widths with potentials xQM and DNR, in keV. Experimental data 



are from [21]. "V. R." means application of the original van Royen-Weisskopf formula, "VRM." 
means the formula taking into account the modifications given by Eq. (|l5|), and "dressed" means 
that a gaussian distribution have been used to dress the quarks. The five last transitions (*) have 
a theoretical width T 



r ee , but the transitions to p p have not been observed experimentally. 







DNR 


X QM 


DNR 


X QM 


DNR 


X QM 




transition 




V. R. 


V. R. 


VRM. 


VRM. 


dressed 


dressed 


Exp. 


p — > e + e~ 




11.360 


0.998 


8.601 


1.315 


5.719 


3.648 


6.75 ± 0.33 


-> p+p- 




11.335 


0.996 


8.582 


1.312 


5.707 


3.640 


6.91 ± 0.42 


u — > e + e~ 




1.217 


0.160 


0.922 


0.187 


0.613 


0.460 


0.60 ± 0.02 


-> p+p- 




1.215 


0.160 


0.920 


0.187 


0.612 


0.459 


<1.52 


(j) — ► e + e~ 




2.984 


0.982 


2.464 


0.964 


1.643 


1.214 


1.30 ± 0.03 


-> 




2.982 


0.981 


2.462 


0.963 


1.642 


1.213 


1.65 ± 0.22 


JMlS)-> 


e + e 


7.381 


11.461 


6.634 


9.740 


5.330 


4.346 


5.16 ± 0.31 




p+p- 


7.381 


11.461 


6.634 


9.740 


5.330 


4.346 


5.11 ± 0.31 


V>(2S)-> 


~e~ 


3.696 


4.603 


3.101 


3.697 


2.074 


1.075 


2.44 ± 0.45 


-» A* 




3.696 


4.603 


3.101 


3.697 


2.074 


1.075 


2.85 ± 1.02 


— > T 




1.436 


1.788 


1.204 


1.436 


0.806 


0.418 


not seen 


^(4040)-> 


e + e~ 


2.681 


2.904 


2.175 


2.273 


1.300 


0.507 


0.73 ± 0.25 


-> 


p+p- 


2.681 


2.904 


2.175 


2.273 


1.300 


0.507 


seen 


-» 


T + T - 


1.768 


1.915 


1.434 


1.500 


0.858 


0.334 


not seen 


V>(4415)-> 


e + e~ 


2.087 


1.979 


1.653 


1.543 


0.900 


0.287 


0.47 ± 0.24 


-» 


p+p- 


2.087 


1.979 


1.653 


1.543 


0.900 


0.287 


not seen 


-> 


T + T - 


1.640 


1.554 


1.298 


1.212 


0.707 


0.225 


not seen 


V>(3770)-> 


e + e~ 








0.023 


0.019 


0.017 


0.009 


0.26 ± 0.05 


-> 


p+p- 








0.023 


0.019 


0.017 


0.009 


not seen 


-» 


T + T - 








0.011 


0.009 


0.008 


0.004 


not seen 


V>(4160)-> 


e + e~ 








0.032 


0.023 


0.021 


0.009 


0.78 ± 0.37 


-» 


p+p- 








0.032 


0.023 


0.021 


0.009 


not seen 


-» 


T + T - 








0.023 


0.016 


0.015 


0.006 


not seen 


T(lS)-> e" 




1.415 


5.052 


1.299 


4.291 


1.287 


2.727 


1.25 ± 0.07 


-» M 


+ M~ 


1.415 


5.052 


1.299 


4.291 


1.287 


2.727 


1.30 ± 0.06 


— ► T 


+r- 


1.404 


5.012 


1.288 


4.257 


1.277 


2.706 


1.40 ± 0.09 


T(2S)-> e" 


^e" 


0.654 


1.496 


0.588 


1.252 


0.583 


0.742 


0.52 ± 0.12 


-» A* 


+ M" 


0.654 


1.496 


0.588 


1.252 


0.583 


0.742 


0.58 ± 0.13 


— > T 




0.650 


1.487 


0.585 


1.244 


0.580 


0.737 


0.75 ± 0.71 


T(3S)-» e 




0.488 


0.922 


0.432 


0.765 


0.428 


0.446 


seen 


-» A* 




0.488 


0.922 


0.432 


0.765 


0.428 


0.446 


0.48 ± 0.08 


— > T 




0.485 


0.917 


0.430 


0.761 


0.426 


0.444 


not seen 


T(4S)-> e" 




0.402 


0.682 


0.353 


0.556 


0.350 


0.326 


0.39 ± 0.17 


-* M 




0.402 


0.682 


0.353 


0.556 


0.350 


0.326 


not seen 


— > T 


+ T- 


0.400 


0.679 


0.352 


0.553 


0.348 


0.324 


not seen 


T(5S)-» e 




0.351 


0.528 


0.306 


0.424 


0.303 


0.252 


0.31 ± 0.08 


-» /i 




0.351 


0.528 


0.306 


0.424 


0.303 


0.252 


not seen 
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0.349 


0.526 


0.305 


0.422 


0.301 


0.251 


not seen 


T(6S)-» 






0.208 


0.432 


0.179 


0.347 


0.161 


0.207 


0.13 ± 0.05 




M' 


V 


0.208 


0.432 


0.179 


0.347 


0.161 


0.207 


not seen 




r" 




0.207 


0.430 


0.178 


0.346 


0.160 


0.206 


not seen 


w(1420) 






0.363 


0.039 


0.233 


0.033 


0.094 


0.041 


seen * 


p(1450) 




e + e~ 


3.063 


0.202 


1.970 


0.192 


0.794 


0.326 


seen * 


w(1650) 




e + e~ 








0.021 


0.012 


0.011 


0.009 


seen * 


0(1680) 




e + e~ 


1.010 


0.280 


0.722 


0.231 


0.278 


0.172 


seen * 


p(1700) 












0.179 


0.094 


0.091 


0.077 


seen * 
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